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Introduction 

We have simulated the performances of a GEM with a large density of multiplication holes. 
The elementary cell is an equilateral triangle whose side is 90/um. We shall assume that this 
pattern extends in the (x,y) plane. At each vertex of the equilateral triangle there is a GEM 
hole with an external radius of 30/um and an internal radius of 20/im. The reference frame 
used in this study has the origin of axis in the center of the GEM hole with the z axis pointing 
to the drift plane. The geometry of the GEM is shown in figure ^ 

In this short note we will describe the simulation of the GEM and, in particular we will 
study the gain and diffusion of the charge for different gas mixtures. This study has been 
performed to finalize the design of the PixiE Imager Detector. 

We have started the simulation by generating single electrons in different positions in the 
(x,y) plane at fixed z-coordinate. This is the most elementary element through which we can 
simulate tracks and the imaging performance of the detector. The process has been followed 
through multiplication in the large fields of the GEM and diffusion of produced electrons 
reaching the readout plane. Due to the cylindrical symmetry of the GEM the first quadrant 
(x > 0,y > 0, see figure ^) has been selected to produce the coordinate (r,(p) of the starting 
electrons at the quota of 40 fim (approximately 15 fim over the top GEM): 

r = 10 -nym (n = 1, 2, 3) <f> = ^ • ^ ~ ^ (k = 1, 2, 3, 4) 

Where 4> is the azimuthal angle. 

At each point we have generated 25 events. The study has been performed for the following 
gas mixtures: 

• 100%CO 2 atm 0.5 ed 1 atm. 

• 20%Ar/80%DME, 50%Ar/50%DME, 80%Ar/20%DME ad 1 atm. 
* http://glastserver.pi. infn.it/pixie/pixie. html 
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Figure 1: A section of the hole of GEM. The red circles indicate the starting points of electrons. 

Gain Study 

We have defined as absolute gain the number of electrons which reach the quota z = —AOfim 
(below the plane the GEM, approximately 15/um below the bottom GEM plane). At this 
quota most of the multiplication processes at the GEM hole are done. 

However not all these electrons drift to the readout plane, some recombine and many stick 
to the lower GEM plane (re-attachment). For this reason we have defined also an effective 
gain as the number of electrons which arrive at the quota z = — 200/im. The electrons 
reaching this quota are considered to be collected by the read out plane. 

It's customary to describe the gain with a Polya: 



Where b is an adjustable parameter and n is the average gain. We have used this formula 
and fit the data with the function: 



where p 2 is a normalization factor and the gain is the Polya mean, G = po • p\. 

The gain distribution is different for different mixtures of gases. In particular the absolute 
gain distribution is often wide with long tails and a description with a single Polya is not 
always satisfactory. Hence, we have described the distribution with the sum of two Polya, 
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of which the first one fits most of events and the second one accounts for the long tails. An 
example is shown in figure [2 We have taken as the average gain of the GEM the mean of the 
first Polya. 
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Figure 2: Distribution of absolute (left) and effective (right) gain, the gas mixture is 50% Argon %50 
DME. The GEM is operated of 480^ with a collection field of AKV/cm. 



Sometimes the distribution shows two clear maxima and a two Polya fit is satisfactory. In 
this case the mean of the two Polya is the average gain of the GEM under analysis. Results 
are shown in the figure 01 0] for a collecting field E t = AKV/ cm and two gas mixtures and in 
the tabled for a collecting field of E t = 5KV/cm. 
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Figure 3: Dependence of effective gain on the GEM voltage. The simulation concerns a gas of 100% 
CO2 and a collection field of AKV/cm. 



The GEM gain increases with the voltage different GEM according to an exponential curve. 

Diffusion Study 

The study of the diffusion of the charge in the collecting region of the detector is important 
for two different issues. 
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AV G em = 560 V E Trasf = hKV/cm 


p = Q.batm 


p = latm 


Gain e ff 


830 ± 100 


70 ± 10 


Gain ass 


7420 ± 400 


2400 ± 900 



Table 1: Average of effective and absolute gain of 100% CO2 at 0.5 and 1 atm and for AVqem of 
560 V and collection field of 5KV/cm. 
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Figure 4: Left: gas mixtures of 20%Ar-80%DME and 50%Ar-50%DME. The plots show the depen- 
dence of the effective gain on the GEM voltage. The collection field is AKV/cm. Right: gas mixture 
80%Ar-20%DME with a large production of secondary avalanches. The figure shows the mean of 
effective gain as a function of the GEM voltage and for a collection field of 4KV/cm. The values of 
gain are the mean of first (gain low) and second (gain hight) Polya which fit the two peaks observed 
in the gain distribution. 



Firstly to establish if the GEM keeps 
memory of the starting point of the elec- 
trons, both in azimuth and radius (with re- 
spect to the center of the hole where the 
avalanche occurs) with a better resolution 
than the granularity of GEM's hole. 

For this, we have studied the position 
of the barycentre of charge arrived on read- 
out plane (barycentre of the avalanche) as 
a function of the position of the starting 
point. 

The second point is the RMS of Gaus- 
sian distribution of charge in the collection 
gap which is related with the spatial reso- 
lution of the detector. 

The average position of the collected 
charge indicates where the multiplication o 
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Figure 5: AR in function of 
at the GEM hole. 
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We have considered 

AR = Rfi n — Rin 

where Ri n = r and Rf% n is the radius of the average charge at the quota z = — 200/xm. 
First of all we have studied the dependence of AR on Ri n 

The figure shows an example of such a dependence. AR is a linear function of Ri n and 
a parameterization: AR = P2Rin + Pi with P2 ~ —1 and Pi « 20u:m is a good fit for 
all simulations (figure EJ). For a ideal GEM AR = in fact this means that the position (in 
radius) of the multiplied charge is the same as the initial electron (the GEM does not disturb 
the image). The results of the fit shows instead that the average collected charge position 
is independent of the starting position (P2 ~ —1) and that the multiplication occurs at the 
radius r 20 — 30«m (the lower external radius). 
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Figure 6: The values of pi (/xm) and p2 in function of percentage of Argon and Neon. 



The last step is the study of the dispersion of the avalanche's charge after it drifted to the 
collection electrodes. 

We have averaged the RMS of the events produced at each point and verified that its value 
is independent of the position of primary electron. Hence we have averaged the RMS of all 
events at all points, to improve the statistics and studied its dependence on the GEM voltage. 
Since, again, we have found no dependence, we have average on all events for a defined gas 
composition. 

The results for 100%CO2 gas at 0.5 atm are: 

• RMS P = 12.5 pm 

• RMS sp = 17.67 urn 

The results for Argon and Neon mixtures are shown in figured the RMS decreases mildly 
with increasing percentage of Argon and Neon. 
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Figure 7: Left we show the projected RMS (fJ,m) averaged on all events for two gas mixtures while 
on the right the spatial RMS (fim). 



to perform reliable simulations. 



6 



